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We theoretically study spin and charge currents through finite-size two-dimensional s-wave superconduc¬ 
tor/uniform ferromagnet/s-wave superconductor (S'/_F/5') junctions with intrinsic spin-orbit interactions (ISOIs) 
using a quasiclassical approach. Considering experimentally realistic parameters, we demonstrate that the com¬ 
bination of spontaneously broken time-reversal symmetry and lack of inversion symmetry can result in sponta¬ 
neously accumulated spin currents at the edges of finite-size two-dimensional magnetic S/F hybrids. Due to 
the spontaneous edge spin accumulation, the corners of the F wire host the maximum spin current density. We 
further reveal that this type edge phenomena are robust and independent of either the actual type of ISOIs or ex¬ 
change field orientation. Moreover, we study spin current-phase relations in these diffusive spin-orbit coupled 
S/F/S junctions. Our results unveil net spin currents, not accompanied by charge supercurrent, that sponta¬ 
neously accumulate at the sample edges through a modulating superconducting phase difference. Finally, we 
discuss possible experimental implementations to observe these edge phenomena. 

PACS numbers: 74.50.+r, 74.45.+C, 74.25.Ha, 74.78.Na 


I. INTRODUCTION 

Spintronics devices operate by spin transport 
mechanisms'"* rather than by utilizing charged carriers, 
as is done typically in conventional electronics devices. 
The use of spin currents can result in higher speeds and 
reduced dissipation^’* while exhibiting weak sensitivity to 
nonmagnetic impurities and temperature."'"* For functional 
spin-based devices, it is necessary to manipulate and generate 
spin-currents in a practical and efficient manner. For this 
reason, many investigations have focused on harnessing 
the spin-orbit interactions'^’*-' (SOIs) present in many ma¬ 
terials, including semiconductors.^""’*"'"*^ The SOI is a 
quantum relativistic phenomenon that can be divided into 
two categories: i) intrinsic (originating from the electronic 
band structure of the material) and ii) extrinsic (originating 
from spin-dependent scattering of impurities)."*"* The in¬ 
trinsic spin-orbit interactions (ISOIs) such as Rashba*^ and 
Dresselhaus*^, are experimentally controllable via tuning a 
gate voltage**’**"^"*. This particular attribute has proliferated 
efforts striving for high-performance spin-based devices, 
including transistors, and new routes in information storage 
and transport. ^’**’**’^*’’^*"^^ Similarly, ferromagnet (F) and 
superconductor (S) heterostructures have received renewed 
interest lately due to the possibility of generating spin polar¬ 
ized triplet supercurrents^*"^^ that can be used for practical 
purposes"^. By considering a ferromagnet with an ISOI, the 
spin orbit interaction can couple with the magnetic exchange 
field, resulting in modified superconducting proximity effects 
and additional venues for new spin phenomena. Indeed, the 
ISOI can induce long-range proximity effects in uniformly 
magnetized S/F structures due to the momentum-dependence 
of the effective exchange field^^. It is therefore of funda¬ 
mental importance not only to find a simple, experimentally 
accessible structure that can support spin currents in F/S sys¬ 
tems, but it is also crucial to determine the spatial behavior of 


the spin currents near the boundaries of the superconducting 
hybrids. 

Many past works are based on the application of external 
electric or magnetic fields. One of the earliest such cases in¬ 
volved the combination of SOIs and an external electric field, 
giving rise to an accumulation of spin currents at the edges""’"*^ 
(the so called spin-Hall effect^®’"***). The spin currents gen¬ 
erally tend to peak near the sample boundaries and vanish 
at the electrode/sample interfaces.These theoretical pre¬ 
dictions were later experimentally observed in semiconduc¬ 
tor samples"*-*. The spin-Hall phenomenon was also exten¬ 
sively studied in superconducting heterostructures where var¬ 
ious types of spin-orbit coupling (SOC) play key roles. ^^’***“*^ 
For example, the out-of-plane component of the spin den¬ 
sity was theoretically investigated*** in a spin-orbit coupled 
S'/W/S'junction [with normal metal {N) interlayer] subject to 
an inhomogeneous external magnetic field. It was found that 
the spin density varies along the transverse direction, leading 
to a longitudinal phase difference between the S electrodes. 
The influence of extrinsic SOIs on the critical supercurrent 
in diffusive magnetic hybrid structures was also studied. *'*’^^ 
In nonmagnetic S/NIS Josephson junctions with SOC sub¬ 
ject to a magnetic field, O-tt transitions may be induced.^® In 
an earlier work*"*, singlet-triplet pair conversion was numeri¬ 
cally investigated using a lattice model in a ballistic half metal 
ferromagnetic Josephson junction with an interfacial Rashba 
SOC. Several optimal configurations have also been theoreti¬ 
cally studied for generating and detecting the predicted long- 
range triplet correlations in experiments.^* 

In this paper, we study the local spin currents in uniformly 
magnetized S/F/S Josephson junctions with spatially uni¬ 
form intrinsic SOIs, avoiding any external electric or magnetic 
fields. We employ a two-dimensional quasiclassical Keldysh- 
Usadel approach that incorporates a generic spin-dependent 
vector potential to study the behavior of the spin current com¬ 
ponents. We consider two types of ISOCs: i) Rashba and 
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ii) Dresselhaus SOC, and the magnetization of the F wire 
can take arbitrary orientations. We hnd that the coupling of 
the F wire’s exchange field and ISOIs leads to edge spin cur¬ 
rents with three nonzero components. The spin current density 
components peak near the edges of the F strip and sharply 
decline when moving away from the edges. Therefore, the 
maximal spin current accumulation takes place near the F 
wire’s corners. This phenomenon can be also observed in ISO 
coupled S'/TV/S'junctions with a single spin active N/S inter¬ 
face as demonstrated in Ref. 75. Our results show that the 
spin current can be switched on or off at the S/F contacts, 
depending on the magnetization direction. The spatially av¬ 
eraged spin current components reveal a 27r-periodicity and 
even-functionality in ip, the phase difference between the S 
terminals. This is in contrast to the charge supercurrent which 
is a 27r-periodic odd function of (p (and consistent with pre¬ 
vious studies®^). Note that for such junctions, the argument 
in the current phase relation for some situations can become 
modified by a simple po shift. 

The simple hybrid structure proposed here relies solely on 
the intrinsic properties of the system, in contrast to other struc¬ 
tures that rely inextricably on external fields to observe the 
edge spin currents^®^®’"^®’^®’®^. Our device consists of a hnite- 
size intrinsic SO coupled F wire (with uniform magnetiza¬ 
tion) sandwiched between two S banks. The spin currents 
then spontaneously accumulate at the sample edges, without 
the application of an external electric or magnetic field to the 
system. We demonstrate that the device allows for the realiza¬ 
tion of spin currents in the absence of charge supercurrent by 
simply modulating p. The edge spin accumulation is a signa¬ 
ture of the spin Hall effect®^, and hence can be experimentally 
measured by optical experiments for instance"^^ (see the dis¬ 
cussions in Sec. III). Also, we discuss the symmetries present 
among the spin current components when varying the magne¬ 
tization orientation with Rashba or Dresselhaus SOC present. 
Moreover, we compare our results with the charge and spin 
currents found in a nonmagnetic diffusive S/N/S Josephson 
junction with Rashba and/or Dresselhaus SOC. We hnd that 
the spin currents vanish in the S/N/S devices, consistent with 
previous works®^’^^, and that the charge current displays a spa¬ 
tially uniform prohle without any transverse component, indi¬ 
cating conservation of charge current. 

The paper is organized as follows. We outline the theoreti¬ 
cal framework used to study hybrid structures with ISOCs in 
Sec. II. In Sec. Ill, the results of diffusive S/F/S Josephson 
junctions are presented for the case of Rashba ISOC. We next 
utilize the symmetries in the spin currents to give a simple 
prescription for hnding the corresponding results for the Dres¬ 
selhaus spin-orbit interaction. We hnally present concluding 
remarks in Sec. IV. 


II. THEORETICAL FORMALISM 


The intrinsic SOI is a consequence of the moving carriers’ 
spin interaction with an atomic potential V{r). Therefore, 
the total Hamiltonian of a moving electron in such an atomic 


potential can be expressed as,^®’^® 

where uiq = O.SlMev, is the free electron mass and and c 
is the velocity of light in vacuum. We define V to represent 
the electron’s momentum vector, and d = {dx,dy,dz). The 
vector of Pauli matrices, t, is given in Appendix A. It has been 
shown that the linearized SOC term can be simply accounted 
for as an effective background held that follows SU(2) gauge 
symmetries."^^’^^’^^’^* Hence, it is sufficient to replace partial 
derivatives, appearing in the quasiclassical formalism, by their 
covariants."^^’^^ Another advantage of the SU(2) approach is 
the convenient dehnition of physical quantities such as spin 
currents.™ 

We start with the Usadel equations that enable us to study 
the charge and spin transport through diffusive F/S systems 
with the ferromagnetic regions having arbitrary magnetization 
patterns h{r) = hy{r), 


d,G{r)[d,G{r)] 


ep 3 -f diag[7/(r),7/'^(r)],G(r) , 

( 2 ) 


'H[r) — h{r) ■ a, 


where p and a denote vectors comprised of 4 x 4 and 2x2 
Pauli matrices (see Appendix A), and D represents the diffu¬ 
sive constant of the ferromagnetic medium. We have denoted 
the quasiparticles’ energy by e which is measured from the 
Fermi surface ep- Throughout this work, we focus on the 
low proximity limit of the diffusive regime^*. In this limit, 
the normal and anomalous components of the Green’s func¬ 
tion can be approximated by, F^°{r) ~ 1 and F(r) ^ 1, 
respectively. Thus, the advanced component of total Green’s 
function, G{r), takes the following form; 


G^ir,e) 


( -1 -F{r,-£)\ 

\F*{r,e) 1 )' 


(3) 


where each entry stands for a 2x2 matrix block. Considering 
the Taylor expansion, the advanced component can be given 
by: 

G^{r,s) = 

/ -1 0 -fnir,-£) -f-{r,-e)\ 

0 -1 -f+ir,-e) -4;(r,-e) 

f*f{r,£) f*-{r,£) 1 0 ■ ^ ^ 

\fUrN) ftiir,£) 0 1 y 

Here we restrict our calculations to the equilibrium situa¬ 
tions where the Retarded and Keldysh blocks of total Green’s 
function are obtained by: G^{r) = —{p3G'^(r)/33}^, and 
G^{r) = tanh(e/cBT/2){G'^(r) — G'^(r)}. Here, kp and 
T denote the Boltzmann constant and system temperature, re¬ 
spectively. 

The Usadel equation, Eq. (2), leads to sixteen coupled com¬ 
plex partial differential equations in the low proximity limit 
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that become highly complicated with the presence of intrin¬ 
sic SOI terms. Unfortunately, the resultant system of cou¬ 
pled differential equations can only be simplified and decou¬ 
pled under very limiting conditions,leading to analyti¬ 
cal results. However, for the systems considered in this pa¬ 
per, numerical methods are the most appropriate, and often 
the only possible routes to investigate the relevant transport 
properties.The differential equations must be supplemented 
by the appropriate boundary conditions to properly capture the 
transport characteristics of S/F/S hybrid structures. We thus 
employ the Kupriyanov-Lukichev boundary conditions at the 
S/F interfaces^' and control the intensity of induced proxim¬ 
ity correlations using the barrier resistance parameter, 

C{G(r)5G(r)} • n = [Gbcs( 0),G(r)], (5) 

where n is a unit vector, directed perpendicular to a given 
interface. The solutions to Eqs. (2) for a bulk, even-frequency 
s-wave superconductor reads. 


rR (Q\-( lcoshi?(£) iCT2e*®sinhi?(£) \ 

Best ) 1(726“*®sinh■(?(£) —lcoshi?(£) y’ ^ 

in which, 

’die) = arctanh(^- ^), 

£ 


is defined in terms of the superconducting gap A. Here the 
macroscopic phase of the bulk superconductor is denoted by 
0 , so that the difference between the macroscopic phases of 
the left and right S electrodes are given by 9i — 9r = For 
more compact expressions in our subsequent calculations, we 
define the following piecewise functions; 


s(£) = 6*® Sinh'!?(£) = 




I - A 2 

c{e) = coshi9{£) = 

I 6 ' 


•\/ 6 ^ — A^ 


0 (£^ - A") - 


le 






where 0(x) stands for the usual step function. It is clear 
that the general boundary conditions given by Eq. (5) do not 
permit current flow through the hard wall boundaries of the 
finite-size two-dimensional SIFIS Josephson junction, shown 
in Fig. 1. 

To study the influence of differing types of ISOI on the 
system transport characteristics, we adopt a spin-dependent 
tensor vector potential A(r) = (^Ax{r), Ay{r),Az{r)^, as 

f0ll0WS;47.73.75, 77,78 


Mr) = l{A:{r)T^+Ayir)r^+AM)r^}, M 

Mr) = \[Al{r)T^+Al[r)Ty (7b) 

Az{r) = \[Al{r)T^+Al{r)Ty+Al{r)T^]. (7c) 



FIG. 1. (Color online) Schematic of a finite-size two-dimensional 
magnetic SIFIS Josephson junction. The superconducting elec¬ 
trodes and rectangular ferromagnetic nano-wire are labelled S and 
F, respectively. We assume that the quasiparticle current experiences 
an intrinsic spin-orbit interaction (ISOI) solely inside the F region. 
The thickness and width of the ferromagnetic strip are labeled dp 
and Wp, respectively. The junction is located in the xy plane and 
the SIF interfaces are along the y axis. The F region has a uni¬ 
form exchange field denoted by h and can take arbitrary orientations 
{K^,K«,h^). 


Using the above vector potential, we define the covariant 
derivatives by; 

d = dl — ieA{r). ( 8 ) 


Accordingly, the brackets seen in the Usadel equation, Eq. (2), 
and the boundary conditions, Eq. (5), (as well as the charge 
and spin currents that shall be discussed below, Eqs. (10) and 
( 11 )) take the following form: 

[5,G(r)] = dG{r) - ie[Alr),G{r)]. (9) 


The spin and charge currents are key quantities that lend in¬ 
sight into the fundamental system transport aspects that pro¬ 
vide valuable and crucial information for nanoscale elements 
in superconducting spintronics devices, as described in the in¬ 
troduction. Under equilibrium conditions, the vector charge 
(J'*) and spin current densities can be expressed by the 
Keldysh block as follows:"'^’^® 


J^(r,v7) = 


^+oo I 

(i£Tr|p3(G(r)[a,G(r)])'^}L (10) 


r^rM = Jo 


/ -\-oo 

d£Tr| psA (G(r) [9, G(r)]) ^ | 

-OO ^ 


( 11 ) 


where J§ = N^eD/A, Jq = hjQ/2e, and Nq is the number 
of states at the Fermi surface. The vector current densities 
determine the local direction and amplitude of the currents as 
a function of coordinates inside the F strip. In other words, 
J(r), provides a spatial map to the currents inside the sys¬ 
tem. We designate 7 = x,y, z for the three components of 
spin current, J®'''. The matrices we use throughout our deriva¬ 
tions are given in Appendix A. To obtain the total Josephson 
charge current flowing through the magnetic strip, an addi¬ 
tional integration over the y direction should be performed on 
Eq. (10) (see Fig. 1). The spin-dependent fields yield lengthly 
and cumbersome expressions, the details of which are not pre¬ 
sented here for clarity. Having now outlined the theoretical 
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approach utilized in this paper, we can now present our find¬ 
ings in the next section. 


III. RESULTS AND DISCUSSIONS 

In our computations below, we consider a uniform and 
coordinate-independent vector potential, A(r), i.e. d-A{r) = 
0, so that the spin vector potential is constant in the entire F 
region. A specific choice for the constant spin vector poten¬ 
tial that results in Rashba (a)^^ and Dresselhaus (/3)*^ types of 
SOC is, 

A% = -Ayy = 2/3, 

Ay = -A^y = 2a, 

Al = Al = 0 , ^ ^ 

Az _ AX _ Ay _ ^ 

By substituting the above set of parameters into Eqs. (7), we 
arrive at. 

Ax = /3t“ - ar^, (13a) 

Ay = ar^ - /3r^. (13b) 

The Rashba SOI^° can be described through spatial inversion 
asymmetries while the Dresselhaus SOI^^ is described by bulk 
inversion asymmetries in the crystal structure.Crystallo¬ 
graphic inversion asymmetries^^ or lack of structural inversion 
symmetries'®’*^’*^’*’ in heterostructures may cause the ISOIs 



FIG. 2. (Color online) Spatial profile of the spin current in 
a uniformly magnetized S/F/S Josephson (see Fig. 1) junction 
without ISOI. The magnetic exchange field is oriented along z, 
h — (0,0, and therefore, solely the 2 component of spin 
current J^^{x,y) is nonvanishing. The junction length and width 
are set to dp = 2.0^s and Wf = 2.0^s, respectively. The 
top panel exhibits the spin current variations along the ^-position 
(the junction length) at four differing locations along the junction 
width, y — 0.25^s, O.S^s, 0.75^3, l.O^s. The bottom panel shows 
JAtx.y) as a function of i/-position along the junction width, at 
. = 0 . 25 ^ 3 , 0 . 5 ^ 3 , 0 . 75^3 10 ^ 3 . 


considered here. For example, strain can induce such inver¬ 
sion asymmetries"**’*’’®’’^* and thus, ISOIs, or the adjoining 
of two differing materials may generate the requisite interfa¬ 
cial sOIs'®'*’’’*’*"*’*’. Nonetheless, there is no straightforward 
method to measure SOIs in a hybrid structure. One possible 
approach would be first principle calculations** in conjunc¬ 
tion with spin transfer torque experiments’*’*®’*’. The intrin¬ 
sic SOIs are often given by the first-order quasiparticle mo¬ 
mentum, which is locked to their spins. This linearized ap¬ 
proach is a simplification to the more generic picture deal¬ 
ing with higher orders of momentum,*'’*’’®’^®* which can be 
observed in e.g., engineered materials.®’’®* We here assume 
that ISOCs can be described by linear terms in the carriers’ 
momentum.**’’® Candidate materials to support spontaneous 
broken time-reversal and broken inversion symmetries include 
electron liquids with ISOIs, which naturally tend to have a 
Stoner-type magnetism at low densities, and a magnetically 
doped topological insulator surface (or by directly coating a 
topological insulator surface with magnetic insulators).**'®® 
Other promising candidates involve the ferromagnetic semi¬ 
conductors (Ga,Mn)As, where both the electronic structure 
and inherent magnetism make these materials well suited for 
experimental studies.’’*’®* Our quasiclassical approach allows 
us to study systems involving nontrivial magnetizations and 
spin vector potentials with arbitrary spatial patterns. We thus 
consider a finite-sized, uniformly magnetized F wire whose 
exchange field can take arbitrary orientations. In order to de¬ 
termine systematically the behaviors of the spin and charge 
currents, we consider three orthogonal magnetization direc¬ 
tions, namely along the x, y, and 2 axes. In addition, we 
incorporate pure Rashba {a 0, (3 = 0) and Dresselhaus 
(/3 7 ^ 0,a = 0) SOCs that allow isolation of their effects 
relative to the physical quantities under study. When finding 
solutions to the Usadel equation, Eq. (2), and the correspond¬ 
ing current densities [Eqs. (10) and (11)], we have added a 
small imaginary part, 6 « O.OlAg, to the quasiparticles’ en¬ 
ergy, 3 —e -F *(5, to enhance stability of the numerical solu¬ 
tions. The imaginary part can be physically viewed as ac¬ 
counting for inelastic scatterings.®* Due to the presence of 
the finite parameter S, we take the modulus of the currents in 
Eqs. (10) and (11). We normalize the quasiparticles’ energy, 
e, and exchange field h by the gap, Ag, at T = 0. Also, all 
lengths are measured in units of the superconducting coher¬ 
ence length ^s- In our computations, we adopt natural units, 
so that ks = fi = 1. 

To begin, we consider for comparison purposes, an 
S/F/S Josephson junction in the absence of SOCs.’*’’® The 
schematic of the S/F/S structure is depicted in Fig. 1. The pa¬ 
rameters C = 4, \h\ = IOAq and dp = 2.0^s, ensure the va¬ 
lidity of low proximity limit considered throughout the paper. 
To have absolute comparisons, we set h = (0,0, h^) and com¬ 
pute the charge and spin currents using Eqs. (10) and (11), re¬ 
spectively. Figure 2 exhibits the spatial map of the spin current 
for Wp = 2.0^s (see Fig. 1). Since the magnetization orien¬ 
tation is fixed along the 2 direction, Jx^{x, y = yo) is the only 
nonvanishing component of spin current for a given fixed loca¬ 
tion 7/0 ■ The top panel of Fig. 2 illustrates the spatial variations 
of y = yo) along the junction length in the x direction 
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FIG. 3. (Color online) Spatial behavior of the three spin current components, {x, y), {x, y), and {x, y) in a uniformly magnetized 

Rashba S/F/S junction. The exchange field of the ferromagnetic strip points along the z direction: h = (0,0, h^) (see Fig. 1). The panels in 
the top row show the spin current components, jp{x, y), as a function of x at four differing locations along the junction width: y=0.25^s, 
O.S^s, 0.75^s, and l.O^s. The bottom row exhibits y) versus y at a;=0.25^s, O.S^s, 0.75^s, and l.O^s. 


at differing positions along the junction width: yo = 0.25^s, 
0.5^S, 0.75^s, and l.O^g. The macroscopic phase difference 
between the S electrodes is set at a representative value, i.e.. 
If = 7 r/ 2 . The bottom panel in Fig. 2 shows J^^{x = xo,y) 
as a function of y, at xq = 0.25^s, 0.5^5,0.75^s, l.O^s. The 
results demonstrate that the spin current is y independent in 
such hybrid junctions, namely J^^{x = xo,y) = const, (we 
also have found Jy^{x,y) = 0). In other words, it is ap¬ 
propriate to view this type of system as an effectively one¬ 
dimensional junction. The variation of J^^{x, y) along the x 
direction is a consequence of spin torque transfer, and hence 
the spin current is not a conserved quantity.The spin 
current is maximal at the SIF interfaces and vanishes at the 
middle of junction, x = 1.0^5 = dp 12. This is contrast to the 
charge supercurrent in the F region, which is conserved, and 
thus has a constant value within the entire F strip (not shown). 


To identify some of the salient features in Fig. 2, we con¬ 
sider now a simplified one-dimensional S/F/S system, which 
permits analytical expressions for the spin current density. 
To this end, we linearize the Usadel equation, and incorpo¬ 
rate the Kupriyanov boundary conditions, where the super¬ 
conducting electrodes have strong scattering impurities. We 
also still assume that the magnetization is oriented along z: 
h = (0, 0,h^). Correspondingly, we define the dimensionless 
quantity, X± = 2i{eFh^)/ep, in which ep is the Thouless en¬ 
ergy, and the dimensionless x coordinate, x = x/dp ^ [O 7 !]■ 
After some straightforward calculations, we obtain the follow¬ 


ing expression for the charge current [Eq. (10)]: 

„ . /■+°“ , 2ztanh(e/cBr/2) 

Jl[x, f) = sin f J de -- 

e)]^^A+cscA_ + A_ cscA+^ + [s*(e)]^ 

^A+cschA_ -F A_cschA+^ j^. 


(14) 


The charge current in this case is seen to exhibit the usual 
sin If odd-functionality in the superconducting phase differ¬ 
ence. Likewise, by substituting the solutions into Eq. (11), we 
arrive at the following expressions for the spin current com¬ 
ponents: 


j:,-{x,f) = 0, (15a) 

J:yix,f) = 0, (15b) 







2ta,nh(ekBT/2) 

C^A+A_ 


|[s*(e)]^A+ cosh2iA_cschA_ cosf + [s*(—e)]^ 

^A+ csc^ A_ (cos A_ -F cos f) sin[A_ (1 — 2x)] — 

A_ csc^ A+(cos A+ -F cost/j) sin[A+(l — 2i)]^ -F 
[s*(e)]^^A+ coth A_cschA_(sinh[A_{l — 2i)] — 
sinh2iA_ cos f) — A_csch^A+(coshA+ -F cos f) 
sinh[A+(l — 2i)]|. (15c) 
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Equations (15a)-(15c) clearly demonstrate that the only non¬ 
vanishing component of spin current is which is consis¬ 
tent with the exchange field aligned along From Eq. (15c), 
it is also evident that is an odd function of the coordi¬ 
nate X relative to the middle of the junction (and thus vanishes 
there), and an even function of the phase difference, ip. These 
features are entirely consistent with the numerical results seen 
in Fig. 2. 

We now incorporate Rashba and Dresselhaus SOCs, while 
keeping the magnetization orientation intact along the z direc¬ 
tion. The ISOIs are confined within the F region and are not 
present within the S electrodes. Through exhaustive numeri¬ 
cal investigations, we have found several symmetries among 
the components of spin current (discussed below) at three par¬ 
ticular directions of the exchange field. Due to the symmetries 
available among the spin current components, we focus here 
on Rashba SOC. We emphasize that similar conclusions can 
be drawn for Dresselhaus SOC through the symmetries de¬ 
scribed below. Figure 3 exhibits the spatial profiles for the 
spin current density components, J^yix,y), and 

J^^{x,y). A square ferromagnetic strip is considered, with 
dp = Wp = 2.0^s, and the superconducting phase difference 
is equal to tp = 7r/2. The Rashba SOC coefficient is set to a 
representative value a = 2.0^s, without loss of generality.^^ 
The top set of panels show jp{x,y — yo) [7 = x,y,z] 
as a function of the x coordinate at yo = 0.25^5, 0.5^s, 
0.75^s, and l.O^s. Whereas the bottom panels represent the 
same quantities, but now as a function of y at Xq = 0.25^s, 
0-5^S. 0.75^s, and l.O^g. As seen in Fig. 1, the junction 
length and width are parallel to the x and y axes, respec¬ 
tively. The components J^^{x,y = yo) and J^y{x,y = yo), 
shown in the top row of Fig. 3, demonstrate that these spin 
current densities vanish at the S/F contacts. This finding is 
consistent with previous works involving nonsuperconducting 
heterojunctions"^^”^^’"^^’®^. The z component, J^^{x,y = yo), 
however exhibits opposite behavior, and is nonzero at the S/F 
contacts due to the exchange field, which is oriented along the 
z axis. Similarly, as seen in Fig. 2, J^^{x,y = yo) is finite 
at the S/F interfaces near the S reservoirs. One of the most 
important features of the results is seen in the top panels of 
Fig. 3, where two peaks in jp{x,y) emerge near the S/F 
contacts. We restrict the spatial profiles to 0 < x < dp/2 
and 0 < y < Wp/2, since the results are symmetrical with 
respect to a; = l.O^s = dp/2 and y = l.O^s = Wp/2, 
so that the maxima of J/S{x, y = yo) occurs near the edges 
of the F wire [at cc = 0, and x = dp]. Turning to the bot¬ 
tom row of panels in Fig. 3, we see that Jp{x = Xq, y) are 
nonzero at the vacuum boundaries, y = 0, and y = Wp. Here 
also the largest values in the spin current density components 
take place near the transverse edges of the F wire {y = 0, 
and y = Wp). The magnitude of the spin current densities at 
X = 0.25^s are generally larger than the other x positions, in 
agreement with the results of J/Six, y = yo) shown in the top 
row of panels. 

We now consider the effects of changing the magnetiza¬ 
tion alignment in the ferromagnet. Thus, Fig. 4 represents 
the same Rashba spin-orbit coupled S/F/S junction as in 
Fig. 3, except the magnetization of the F wire is now ori¬ 


ented along the y axis. This specific direction of h leads to 
Jx'^(.x,y = yo) = J^^[x,y = yo) = 0 at the S/F inter¬ 
faces and the spin current densities peak near the edges of F 
wire. The spin current density J/y{x,y = yo) however is 
nonzero at the S/F contacts similarly to J/.^{x,y = yo) in 
the configuration where the magnetization points along the z 
direction (Fig. 3). As mentioned earlier, this nonvanishing be¬ 
havior is directly related to the exchange field direction which 
lies now parallel to the y axis. Examining Jp{x = Xo,y) 
in the bottom row panels in Fig. 4, the maximal values of 
J/S{x = Xo,y) take place near the vacuum boundaries, i.e. 
y = 0, and y = Wp. Our investigations demonstrate sim¬ 
ilar qualitative trends for the components of J/S{x = xo,y) 
when the magnetization resides along the x axis. Note that the 
transverse components of the spin currents are nonzero inside 
the ISO coupled ferromagnetic wire, i.e., Jp{x, y) 7 ^ 0, and 
vanish at the vacuum boundaries (y = 0, Wp). We mainly fo¬ 
cus here on the jp{x, y), since the longitudinal components 
contain the relevant information needed to describe and under¬ 
stand the accumulation of spin current densities at the edges 
of the structures. 

Considering now the previous characterization of the spin 
current components in systems with ISOCs, we schematically 
summarize the spatial maps in Fig. 5 for J^'^(x,y). The 
largest amplitudes of J^'^{x,y) reside near the edges of the 
F strip, i.e. x = 0,dp and y = 0, Wp. We have qualitatively 
marked these regions by light yellow “ribbons”. Therefore, 
the overlap of maximal amplitudes take place near the cor¬ 
ner regions of the F strip. We have marked these areas by 
dashed curves and with a deeper yellow color. The spatial 
profiles found here are qualitatively similar to the existence of 
edge spin currents found in nonsuperconducting heterojunc¬ 
tions with ISOIs,except with one crucial difference: here 
spin accumulation at the edges arises in the absence of an ex¬ 
ternal field. As mentioned in other works®^, the spin accu¬ 
mulation is a signature of the spin Hall effect. Therefore, the 
predicted spin accumulation in this paper might be measurable 
through optical experiments"^^, such as through Kerr rotation 
microscopywhere spatial profiles of the spin polarizations 
near the edges can be imaged. An alternate experimental pro¬ 
posal involves multi terminal devices"^"^’"^^. When transverse 
leads are attached to the lateral edges of a two-dimensional 
S/F/S junction (borders at y = 0, y = Wn in Fig. 1(a)), the 
spin accumulations at the F wire’s edges inject spin currents 
into the leads."^"^’"^^ The transversely injected spin currents into 
the lateral leads in turn may induce a voltage drop between 
the additional leads."^"^’"^^ 

The Josephson effect is a significant example of a macro¬ 
scopic quantum phenomenon and is of fundamental impor¬ 
tance in determining the properties of dissipationless coherent 
transport. Thus, the behavior of the spin currents upon vary¬ 
ing the macroscopic phase difference is crucial to experiments 
and applications utilizing spin-Hall effects and spin transport. 
In Fig. 6 , we therefore study the spin current components as 
a function of the macroscopic phase difference, (p, between 
the S banks. We consider the parameter set used in Fig. 4, 
including h = (0, F',0). In the top set of panels, the spatial 
variations of J/P'{x,y) are plotted at tp = 0, 0.27r, 0 . 47 r, and 


7 



FIG. 4. (Color online) Spatial profiles of the spin current components; y), y), and y) in an S/F/S system. The Rashba 

ferromagnetic wire’s width and length are equal ioWp = dp = 2.0^s. The exchange field of the ferromagnetic strip is fixed along the y 
direction, h = (0, ,0). Top row. spatial behavior of y) along the junction length, x, at t/ = 0.25^s, 0.5^s, 0.75^s, l.O^s. Bottom 

row. spatial variations of jp{x, y) along the junction width in the y direction at x = 0.25^s, O.S^s, 0.75^s, l.O^s. 


I.Ott. We have also chosen a representative position along the 
junction width, corresponding to y = l.O^s = Wp!% which 
simplifies the analysis while maintaining the generality of the 
discussion. Although Jp{x, y) has a minimum at j/ = l.O^s, 
it exhibits the same trends as a function of (f compared to the 
other positions inside the F wire. By increasing the super¬ 
conducting phase difference from 0 to it, the amplitudes of 
the spin current components decrease overall. In the bottom 
set of panels of Fig. 4, we illustrate [x,y))^ as a func¬ 
tion of ip, a.t y = 0.25^s, O.S^s, 0.75^s, and l.O^g. Here, 
we denote the spatial average over the x coordinate from 0 
to dp by {...)x. In order to better visualize the averaged pro¬ 
files, we use logarithmic scales in the vertical axes of the bot- 



dp 


FIG. 5. (Color online) Qualitative illustration of the edge spin cur¬ 
rent densities in a Rashba or Dresselhaus spin-orbit coupled S/F/S 
junction. The light yellow ribbons display edge regions with maxi¬ 
mum spin current densities. The induced spin currents can be con¬ 
sidered as a response of the intrinsic spin-orbit coupled system to 
the presence of an exchange field (the combination of spontaneously 
broken time-reversal symmetries and the lack of inversion symme¬ 
tries). As shown in Fig. 1, the exchange field of the ferromagnetic 
wire is uniform and can take arbitrary directions. The regions that 
carry maximal accumulation of spin currents are qualitatively shown 
by the semicircular regions. 


tom row of panels. As seen, the three components of spin 
current {jp{x,y))x are even-functions of ip, with a period 
of 27r, namely -\- ip) = J‘^^{—ip), n G Z. This is 

contrary to the charge supercurrent which is an odd-function 
of p, i.e., J'^(2n7r + p) = —J‘^{—p) regardless of a finite 
phase-shift These findings are entirely consistent 

with previous studies of S/F/S Josephson junctions with in¬ 
homogeneous magnetization patterns®^. We here remark that 
an additional phase-shift pq may appear in such junctions due 
to the coupling of exchange field and ISOIs.®®’^^’^** Nonethe¬ 
less, the explicit current-phase relations simply undergo a shift 
in According to the current-phase relations, the charge 

supercurrent vanishes at certain p that is quite different than 
the behavior of the spin current components which clearly 
show nonzero values at the same p. Therefore, these differ¬ 
ences in charge and spin currents allows for an examination 
of edge spin currents without any net charge current in an ISO 
coupled F wire sandwiched between two S banks. 

We are now in a position to discuss symmetries that may 
arise among the spin current density components for differ¬ 
ing magnetization orientations in systems with either Rashba 
or Dresselhaus SOCs. Our investigations have found that 
the out-of-plane spin current, J^^{x,y), remains unchanged 
upon exchanging the Rashba and Dresselhaus SOCs, regard¬ 
less of the magnetization orientation. This follows from the 
form of the spin vector potential discussed at the beginning 
of this section. However, this picture changes for the in-plane 
js{x,y} y'^ components. The x and y components of the 
spin current become interchanged when transforming from 
one type of spin-orbit interaction to another. Precisely speak¬ 
ing, by going from Rashba to Dresselhaus spin-orbit coupling. 
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FIG. 6. (Color online) Spin current components for differing values of the superconducting phase difference ip'ma. Rashba spin-orbit coupled 
SI FIS Josephson junction. The top row of panels shows the spatial variations of y) as a function of x along the junction length at 

ip = 0,0.27r, O.dyr, and l.Ovr. The location along the junction width is fixed at the middle of the junction, y = Wrl'i ~ l.O^s. The bottom 
row of panels represents the spatially averaged spin current components over the junction length (denoted by ( {x,y))x) ^P- The average 

is performed along four positions: t/=0.25^s, O.S^s, 0.75^s, and l.O^s. The ferromagnetic wire is a square strip with d_F = W_f = 2.05s, 
and exchange field h — (0, 0). 


one simply needs to exchange indices x and y in the compo¬ 
nents of both the exchange field and the spin current. Oth¬ 
erwise, everything stays the same. By making use of the 
simple transformation rules described, one can easily deduce 
the results of Dresselhaus spin-orbit coupled systems from the 
plots presented for Rashba spin-orbit coupled SIFIS systems 
shown in Figs. 3, 4, and 6. 

To conclude this section, we briefly discuss the importance 
of having a magnetic element in the Josephson junction for the 
effect of spin current edge accumulation to take place sponta¬ 
neously. We thus take the limiting case of ft. = 0 in our pre¬ 
vious calculations above involving SIFIS junctions. Using 
otherwise the same geometrical and material parameters, this 
case was found to produce no spin current, = 0, 

in the presence of Rashba (a ^ 0, f3 = 0) and/or Dressel¬ 
haus (/3 7 ^ 0, a = 0) SOIs. These findings are consistent with 
previous works,where several simplifying approximations 
were employed for Rashba-based SINIS systems. Examin¬ 
ing also the charge supercurrent, y), for both the Rashba 
and Dresselhaus interactions, we observed a uniform spatial 
map for the charge current density for all (p, with Jj: = const., 
and Jy = 0. In other words, the spin-dependent fields cannot 
induce transverse charge supercurrents in a diffusive SINIS 
junction. This is in stark contrast to its ballistic SINIS coun¬ 
terpart, where a transverse charge supercurrent (that is, equiv¬ 
alent to a supercurrent flowing along the y direction in our 
configuration depicted in Fig. 1) was theoretically predicted 
due to the presence of intrinsic SOIs^^. 


IV. CONCLUSIONS 

We have theoretically studied the behavior of spin and 
charge currents in a finite-size two-dimensional SIFIS 
Josephson junction with intrinsic spin-orbit couplings. We 
utilized a two-dimensional Keldysh-Usadel quasiclassical ap¬ 
proach that incorporates a generic spin-dependent vector po¬ 
tential. Our results demonstrate that the combination of a uni¬ 
form magnetization and ISOIs drives the spin currents which 
spontaneously accumulate at the F wire’s edges. The comers 
of the F wire were shown to host the maximum density of spin 
currents. (As demonstrated in Ref. 75, similar edge phenom¬ 
ena can be found in finite-size two-dimensional intrinsically 
spin orbit coupled SINIS junctions with a single spin active 
interface. Additionally, it was shown that maximum singlet- 
triplet conversions take place at the corners of N wire nearest 
the spin active interfaces'^.) Our investigations show that the 
spontaneous edge accumulation of the spin currents are ro¬ 
bust and can exist at all magnetization orientations, indepen¬ 
dent of the actual type of ISOIs. Our investigations have also 
found several symmetries among the spin current components 
upon varying magnetization orientations coupled to a Rashba 
or Dresselhaus SOI. By varying the superconducting phase 
difference, (p, between the S banks, we determined the spin 
and charge currents as a function of phase difference. We have 
found that net spin currents therefore emerge and accumulate 
spontaneously at the edges, in the absence of charge flow, 
when properly modulating p in finite-size two-dimensional 
intrinsically spin-orbit coupled SIFIS hybrid structures. This 
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work can be viewed as complementary to previous studies 
involving edge spin currents in non-superconducting spin- 
orbit coupled structures where externally imposed fields were 
required'^^^^’"^^’^®’®^ . We have shown that remarkably, edge 
spin currents can be spontaneously driven by the coupling 
of intrinsic properties of a system, i.e. spontaneously bro¬ 
ken time-reversal and the lack of inversion symmetries in the 
absence of any externally imposed field. 





We also introduced the 4x4 matrices p = {pi,P 2 , Ps)' 
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Following Ref. 68 , we define r'*', , and po as follows; 


Appendix A: Pauli Matrices 





0 \ _ 
(jl* j ’ Po 




In Sec. II we introduced the Pauli matrices in the spin space 

and denoted them by a = (cr^, a^, cr^), f = (r^, t^), and to unify our notation throughout the paper 7 stands for x, y, z. 
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